Interacting two helical edge modes in quantum spin Hall systems 

Yukio Tanaka^, and Naoto Nagaosa^'^ 
^Department of Applied Physics, Nagoya University, Nagoya, 464-8603, Japan 
^ Department of Applied Physics, University of Tokyo, Tokyo 113-8656, Japan 
^ Cross Correlated Materials Research Croup (CMRC), 
ASI, RIKEN, WAKO 351-0198, Japan 
(Dated: April 22, 2009) 

Abstract 

We study theoretically the two interacting one-dimensional helical modes at the edges of the 
quantum spin Hall systems. A new type of inter-edge correlated liquid (lECL) without the spin 
gap is found. This liquid shows the diverging density wave (DW) and superconductivity (SC) 
correlations much stronger than those of the spinfull electrons. Possible experimental observations 
are also discussed. 
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Quantum spin Hall systeni 
usual band insulator l|, 2, [s], 



'QSHS) is a new state of matter topologically distinct from the 
5|. It is protected by the band gap induced by the relativistic 



picture for QSHS is the two copies of quantum Hall system of up and down pseudospins 
with the opposite chiralities, and hence the edge modes with the opposite direction of the 
propagation for different pseudospins along the boundary of the sample are expected 3]- 
Time-reversal symmetry and Kramer's theorem guarantee the stability of these helical edge 
modes, i.e., prohibiting the backward scattering between the time-reversal pair states and 

n Q 

protecting the crossing of the energy dispersions [8|, |9(. These modes, called the helical 
edge modes, are experimentally observed through the quantized charge conductance 2e^ /h 
in the quantum well of HgTe/(Cd,Hg)Te [4]. These helical edge modes on one edge are 
the "half of the spinfull one-dimensional electrons, because only one chirality is allowed for 
each pseudospin. Therefore, several nontrivial features are expected for this one-dimensional 
system such as the half-e charge near the ferromagnetic domain wall and the robustness 
against the umklapp scattering [8|, 19|]. Several other theoretical aspects have been studied 



very recently 
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In this paper, we study theoretically a new state created by the fusion/recombination 
of two of the helical edge modes, i.e., the two half spinfull electrons, which is different 
from the original one-dimensional electrons. This fusion can be realized by the interaction 
between the edges of the opposite sides of the sample, or one can even design the interacting 
helical edge modes with the opposite helicity (Fig. 1(a)) or same helicity (Fig. 1(b)). 
For a moment, we consider the anti-parallel case (Fig. 1(a)), and the translation to the 
parallel case (Fig. 1(b)) is straightforward and will be discussed later. We employ the 
bosonization method. There are four fields involved in this problem, i.e., right moving up- 
pseudospin electron density pi_R|, and left moving down-pseudospin electron density Pm-, 
for the edge 1, and p2Ri, P2LT, for the edge 2. These densities are related to the phase 
fields of electrons ( which are defined below Eq.(4) ) ^i-f, 021, (t>2], respectively, as 
Pii?T — (Vv^)^a;0iT ^tc. Q, Q, Q (see Fig. 1). This means the spatial variation of the 
phase corresponds to the accumulation/depletion of the electrons. The novel states induced 
by the interaction is realized when some combination of these phases becomes rigid, leading 
to the collective motion of the electrons. We have found that the different combination 
of 0's gets rigid in the two interacting helical edge modes from those realized in the usual 
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(a) Anti-parallel 



QSHS 



(b)Parallel 




FIG. 1: Schematic illustrations of two helical edge modes with (a)anti-parallel case and (b) parallel 
case. 

spinfuU electrons. Namely, the helical edge modes offer an unique opportunity to realize 
a new state by the recombination of the "fraction" of electrons. To reach this conclusion, 
we have classified all the possible forms of the interaction, and found that the inter-edge 
pseudospin flip backward scattering plays an important role in the incommensurate case, 
while also the inter-edge umklapp scattering can be relevant in the commensurate case. It is 
found that an inter-edge correlated liquid (lECL) analogous to bipolaron liquid is realized in 
the incommensurate case, showing the density wave and superconducting correlations both 
of which are more strongly diverging compared with the one- dimensional spinfull interacting 
electrons. Furthermore, there is no spin gap in this liquid state, in sharp contrast to the 
Luther-Emery liquid or the spin gap state in ladder system. Physical properties of this novel 
liquid are described proposing possible experiments. 
The electron-electron interaction is originally given by 

i^int = \jdrj dr'ijl{r)ijl{r')V{r - r')^.'(r')^.(r) (1) 

where r, r' are the three-dimensional coordinates, V{r — r') is the interaction potential, and 
the field operators ipKr), ^/'o-(r) describe the creation and annihilation of the electron with 
spin a at the position r, respectively. Now we restrict to the helical edge modes only, which 
means 

Mr)= Yl 5^ C«a^aa(x)^a(y,2;). (2) 
a=t,ia=l,2 

Here, the index a specifies the edge, a is the spin of the electron, while a is the pseudospin 
in the presence of the spin-orbit coupling, and Caa is the coefficient of the transformation. 
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^aa(a^) represents the field operator of tlie one-dimensional electrons along the edge channel, 
while v^ad/, z) is the wave function of the helical edge modes perpendicular to the edge a. 
We put this expression into Eq.([T]), and neglect the overlap of the wave functions for the 
helical edge modes at different edges, which means that the different edges can not appear 
for the same spatial point r or r'. Also the pseudspin a determines the direction of the 
propagation for each edge. 

When we consider the case where the direct overlap of the wave functions between the 
two edges can be neglected, the form of the interaction is restricted. The Hamiltonian 
H = Hk + -f^int is the sum of the kinetic energy Hk and the electron-electron interaction 
Hint as given by 

Hk = -I j dx[^[^d,m^^ - ^{^d^^n 

+ ^\^d.^2i - (3) 

and 

+ '^l^{x)'^l^{x + a)^2ii{x + a)^2^{x)]6 + h.c.} 
+ g'^ J dx[^l^^\^^ii^2^5 + h.c.] 

+ gsf j dx[^^^l^^^^2i + ^l^^^'^^^2r] (4) 

by right movers (left-movers) and "$21 {^11 and 'I'2t) with 6 = expl—iAkpx). gf, g'j, 
du, g'u, and gsf denote the interaction constant of intra-edge forward, inter-edge forward, 
intra-edge umklapp, inter-edge umklapp, and inter-edge pseudospin-fiip backward scatter- 
ing, respectively. The corresponding Hamiltonian with parallel case can be derived by 
substituting and \l/2t for ^2^ and "^21- 

Fermion operators are expressed by = ?7i| exp(— iA/47r0i|)/v^27r, ^'i^ = 

r/i| exp(iV4vr0u)/V2vr, ^2T = V2^ exp{i^/A^T(f)2^')/^/2^T, and \l'2i = r72i exp(— jVl7r02i)/V2vr, 
with boson fields 02T and (f)2i [15]. rj^, r/i|, 7721 and 772; are Klein factors. 

Now we consider the bosonization of the Hamiltonians Eq.Q and Eq.(jl]) [l^. By intro- 
ducing cpi = 0i| + (p2 = (f)2i + 4>2h ^1 = 01? ~ 0U5 s-nd 62 = (j)2i - the Hamiltoniau 



H can be given by 



^9f_ f ^^^dipidip2 061862^ 



2tt J dx dx dx dx 



27r2 

I 9'u 
27r2 



j dx{cos[A{^ipi — kpx)] + cos[4(y7rv^2 — kpx)]} 
( J dx cos[-\/47r(v9i + (^2) — 4:kpx] 



-^C j dxcos[VA^{ipi- ip2)] (5) 

with C = ?72j,?7it^ij,^72T- III order to diagonalize the inter-edge forward scattering terms {g'f), 
we introduce the symmetric and antisymmetric combinations by (p'^ = {(fi + (P2) / V2, ip'^ = 
{(pi - f2)/V2, e'^ = {61 + 62) /V2, and 6'^ = {Oi - 62) /V2. Then, the above Hamiltonian 
Eq. ([5]) can be transformed into H = Hq + + H^' + Hgf with each of the term being given 
by 



Hu = ^ dx cos{VS7T(p's — 4:kFx) cos^VSnip'^), (7) 

TT^ J 

H^, = ^( I dxcoslVs^if'^-Akpx], (8) 



27C 

Hsf = -0C / dxcosiVSn^',). (9) 



The dimensionless parameters Ki and K2 in Eq. are given by Ki = \ / "^^^^j^-^p^^ and 



K2 = y all+igj-g^) ^i*^ ^1 = ^ + (S-/ + 5'/)/(2vr) and = + (5-/ - 5(})/(27r). Note that 
the difference between Ki and K2 is due to the presence of inter-edge forward scattering 
Now, we consider the relevance/irrelevance of the nonhnear terms H^, Hu' and Hgf. These 
terms are relevant for Ki + K2 < 1-, Ki < 1 and K2 < \, respectively. For the single edge 
case, Ki = K2 = K, becomes relevant for K < 1/2 8|, which requires strong repulsive 
force. The inter-edge interactions Hu' and Hgf can be relevant more easily. The terms if„ 
and Hu' survive only in the commensurate case; Akp = G with a reciprocal lattice vector G. 
Hsf, on the other hand, exists even for the incommensurate case. 
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Let us first consider the incommensurate case, which is the most general case because 
the crossing energy of the helical edge modes are not necessarily at the center of the gap 
without the particle-hole symmetry, and also the Fermi energy is often determined by the 
details of the sample, e.g., in-gap impurity levels and does not give the commensurate filling 
of the edge channels. In this case, only the inter-edge pseudospin-fiip backward scattering 
Qsf is relevant for K2 < 1, which is realized when the intra-edge repulsive forward scattering 
gf{> 0) is stronger than the inter-edge repulsion g'f{> 0) for example. This relevant Qs/ term 
fixes to open the gap for its fluctuation. The remnant massless modes are the conjugate 
pair (fi'g and 6'^, which correspond to the correlated motion of the charge /spin between the 
two edges analogous to the bipolaron formation in spinfull ladder system 16|, ll7| . 

We first focus on the generic phase diagram. It is determined by comparing the critical 
exponents of diverging correlation function of the following four order parameters, i.e., 
(i) Intra-edge Density wave (DW), ^1^^'u±^^^^'2t; (ii)Inter-edge DW, ^l^^'2T±'S'u^2i; (iii) 
Intra-edge Superconductivity (SC) \E'i|^'ij±\l/2j^'2T) and (iv) Inter-edge SC, \I'i|^'2t='=^ii^2i- 
Using the standard bosonization method, the correlation function of each of the above order 
parameters behaves as 

< O{x,t)O{0,0) >~ [max(| x \,Vi=i,2 \ r 1)]"^ 

in the presence of gapless modes. Figure 2 shows the phase diagram in the plane of gj and 
g'jr. The phase boundaries, which are the straight lines gf = g'j, gf = —g'j- correspond to 
K2 = 1, Ki = 1, respectively. Namely, each region is characterized by (A) Ki < 1, K2 < 1, 
(B) Ki > 1, K2 < 1, (C) Ki > 1, K2 > I, and (D) Ki < 1, K2 > 1, respectively. Table 
I and Table II explain the behaviors of the order parameter correlation functions for the 
incommensurate and commensurate cases, respectively. For incommensurate case (Table I), 
Hsf is relevant and lECL is realized in regions A and B. In this case, only the (i) Intra- 
edge DW or (iv) Inter-edge SC shows the critical power-law behavior, while the others are 
exponentially decaying (" — " symbol in Table I). The most dominant order in each region is 
indicated by the bold letters, i.e., (i)Intra-edge DW in region A, (ii) Inter-edge DW in region 
D, (iii) Intra-edge SC in region C, and (iv) Inter-edge SC in region B. 



We further study the physical properties of the lECL other than the order parameters. 
First, we consider the stability of this state against the inter-edge tunneling which has been 
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FIG. 2: Phase diagram of two helical edge modes. Intra-edge Density wave (DW), Inter-edge 
Superconductivity (SC), Intra-edge SC, and Inter-edge DW, are most dominant in the region A, 
B, C, and D, respectively. 



Region 


A 


B 


C 


D 


state 


lECL 


TL 


Intra 
DW 




Ki 


K1 + K2 




Inter 
DW 




Ki + I/K2 


Ki + I/K2 


Intra 
SC 




I/K1 + I/K2 


I/K1 + I/K2 


Inter SC 




1/i^i 


l/Ki + K2 


l/Ki + K2 



TABLE I: The relation between the most dominant order parameters and the exponents of cor- 
relation function for incommensurate case. " — " denotes that the correlation function has an 
exponential dependence. The most dominant order in each region is indicated by the bold letters. 

neglected so far. The Hamiltonian describing this process is given by 

Ht = -lj dx{t cos[V2^{v', + 9'^)] - t'cos[v^(^; + e',)]} (10) 

where t and t' are transfer integrals without and with pseudospin-flip scattering, respectively. 
It is evident that this perturbation is irrelevant since ip'^ is fixed in lECL. 

Next, we consider the magnetic properties of lECL. In this state, the spin-spin correlation 
function has two contributions, < S{x = 0,t)S{x = 0,0) >~ + {^)^^. The Fourier 
transform with respect to the imaginary time r thus leads to the spin susceptibility Xs ~ 
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rp _^rpKi-i ^ function of temperature T 16|. Then the resulting temperature dependence 
of NMR relaxation rate Ti is given by ^ ~ T + T^^^^. It should be remarked that 1/Ti 
has a divergent T dependence as T ^ for Ki < 1, where Intra-edge DW state is dominant 
realized in region A. This anomalous temperature dependence of Ti 181] can be detected by 
experiments of NMR. 

The Josephson coupling through lECL shows also a novel property. Here, we consider the 
lECL with length d is sandwiched by two superconductors Si and S2 [l9]. The Josephson 
current Ij is given by Ij = —2edF{x)/dx with phase difference x- The phase depen- 
dent part of the free energy is proportional to F{x) oc Re[tl{t*2fe-'^ dTU{d,T)] with 
n(d,r) =< O(0,0)O(rf,r) >, 0(0,0) = ^n^'at ± ^u*2i and (3 = 1/T. ^1(^2) is the trans- 
mission amplitude at IECL/Si(S2) interface. Following the discussion by Fazio et al. jigj ]. 
the resulting Ij at zero temperature is proportional to d^^/^^^^ for Ki < 1. This is quite 
unusual that Josephson current is hugely enhanced as compared to conventional spinfuU 
Luttinger chain where Ij is proportional to d~^^^ . 



Region 


A 


B 


C 


D 


state 


DW 


lECL 


TL 


lECL 


Intra 
DW 


Order 


Ki 


K1+K2 


K2 


Inter 
DW 






Ki + l/i^2 




Intra SC 






17X1+ 1/^2 




Inter SC 




1/Ki 


l/Ki + K2 





TABLE II: Similar Table to Table I for the commensurate case. 



Now we turn to the commensurate case (Table II). The term H^' [Eq. ([8])] fixes the phase 
(p'g for Ki > 1 which is realized in regions A and D. The Intra DW phase is stabilized in 
the region A and order parameter has a nonzero value, i.e., long-range ordered state at zero 
temperature. The Inter DW is also stabilized and the resulting u becomes I/K2 in the 
region D. On the other hand, the exponents z/'s in the regions B and C do not change from 
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the incommensurate case, since the umklapp terms is irrelevant there. Note that the most 
dominant order parameter in each region does not change from the incommensurate case. 

In the parallel helical edge modes (Fig. 1(b)), we can obtain essentially the same results 
as the anti-parallel case discussed thus far. The only difference is that the corresponding 
operators for (ii) Inter-edge charge density wave (CDW)/ spin density wave (SDW) and 
(iv) Inter-edge SC are given as \l/|-|.'i'2j, ± ^'J|\l/2t and \l/i|^'2i ± ^i|\E'2T- This means that 
the opposite pseudospin electron operators appear in the parallel case instead of the same 
pseudospin ones in the anti-parallel case. 

Finally, we compare the present lECL with the one-dimensional liquids already known. 
The single edge mode case is reproduced by setting g'j = 0, g'^ = and gsf = 0. Then 
Ki = K2 = K is satisfied. The intra-edge umklapp scattering g^ becomes relevant for 
K < 1/2 [8]. The exponent of correlation function u of Intra-edge SC and Intra-edge DW 
are given hj 2/K and 2K, which shows weaker correlation compared with those in Table 
I and II. For spinfuU one-dimensional TL liquid, the u of Intra-edge SC and Intra-edge 
CDW/SDW are given by 1/Kp + K^{l/K„) and Kp + K„{l/K„) [l6|. As compared to this 
spinfull one- dimensional TL liquid, the density wave (SDW/CDW) and superconducting 
correlations both of which are more strongly diverging. The spin gap phase is realized in 
Luther- Emery liquid [2^. It is known that the v of the intra-edge SC and the intra-edge 
CDW are given by 1/i^' and respectively. Although these exponents are similar to those 
in lECL, lECL is very much different since there is no spin gap. This difference comes from 
the fact that the combination + 02t ~ 02; — 0ii is fixed in Luther- Emery liquid while a 
different one (pi^ — 02t ~ 02i + <Pii is fixed in lECL. The present lECL is also different from 
spinfull two-leg ladder, where the v of the inter-ladder d-wave pairing can become 1 / (2i^ ) 



2l| . However, there opens the spin gap also in this state, in sharp contrast to the present 



lECL. 

In this letter, we have proposed an Inter-edge correlated liquid (lECL) for the two helical 
edge modes of quantum spin Hall systems. lECL is realized most probably in the incommen- 
surate case, showing the density wave and superconducting correlations both of which are 
more strongly diverging compared with the one-dimensional spinfull interacting electrons. 
Furthermore, there is no spin gap in this liquid state, in sharp contrast to the Luther-Emery 
liquid or the spin gap state in two-leg ladder system. This state can be clearly identified 
experimentally by NMR and Josephson junctions. 
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